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Abstract: The propagation of UHECRs is affected by the intergalactic magnetic fields that were pro- 
duced during the course of the large scale structure formation of the universe. We adopt a novel model 
where the large scale extragalactic magnetic fields (EGMF) are estimated from local dynamic properties 
of the gas flows in hydrodynamic simulations of a concordance ACDM universe. With the model mag- 
netic fields, we calculate the deflection angle, time delay and energy spectrum of protons with E > 10^^ 
eV that are injected at cosmological sources and then travel through the large scale structure of the uni- 
verse, losing the energy due to interactions with the cosmic background radiation. Implications of this 
study on the origin of UHECRs are discussed. 
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Introduction 

The enigma of ultrahigh energy cosmic rays (UHE- 
CRs) with energy £; > 5 EeV (1 EeV = lO^^ eV) 
has received considerable attention, both observa- 
tionally and theoretically, over the past decades 
[5, 7]. Although the composition studies indi- 
cate that most of them are likely to be protons 
[1], the nature and origin of these CRs remain yet 
to be understood. At these high energies, UHE- 
CRs are not confined within the galactic magnetic 
plane and therefore their origin should be extra- 
galactic. In particular, the isotropic distribution 
of arrival directions suggests that large number 
of UHECR sources might be distributed at suffi- 
ciently large distances [7, 4]. UHECRs from cos- 
mological sources get significantly attenuated dur- 
ing their propagation due to interaction with the 
cosmic microwave background, resulting in the so- 
called Greisen-Zatsepin-Kuzmin (GZK) cutoff in 
the energy spectrum [3]. Interestingly, AGASA ex- 
periment failed to identify this cutoff [7], whereas 
HiRes recorded it [9]. Detection of super-GZK 
events with isotropic arrival directions is enig- 
matic, because powerful astronomical objects are 
very rare in the local universe. New results from 



Pierre Auger experiment reported in this confer- 
ence may shed light on this puzzle. 

Charged UHECRs are expected to be deflected by 
large scale extragalactic magnetic fields (EGMF), 
so their arrival directions can deviate from source 
directions and arrival time is delayed compared to 
the light travel time [2, 6]. At present, the main 
challenge is to estimate the EGMF, since observ- 
ing it is still a difficult task. From this standpoint, 
we employ a novel model for the EGMF in which 
the strength of the intergalactic magnetic fields is 
estimated from local turbulence in a cosmological 
hydrodynamic simulation. The trajectories of UHE 
protons from active galactic nuclei (AGNs), prop- 
agating through the model EGMF, are calculated, 
including all relevant energy loss processes. We 
then study the angular deflection, time delay and 
energy spectrum of UHE protons that arrive in the 
regions similar to the Local Group. 

Model 

Extragalactic Magnetic Fields 

Concordance ACDM cosmological simulations are 
carried out along with passively evolved magnetic 
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Figure 1 : (a) Volume fraction in the gas density-EGMF strength plane with our model EGMF at z = 0. (b) 
Volume fraction, df /d\og{B), (soUd Une) and its cumulative distribution, /(> B), (dotted) as a function 
of EGMF strength. 



fields in a cubic box of comoving size lOO/i^^ 
Mpc, using 512'^ grid zones with the following pa- 
rameters: D.BM = 0.043, VIdm 0.227, and 
= 0.73, h = i/o/(100 km/s/Mpc) = 0.7, and 
as = 0.8. Six different initial realizations are 
used. Only the directional information of the pas- 
sive fields is adopted. Assuming that the inter- 
galactic magnetic fields result from turbulent mo- 
tion of the intergalactic gas, the strength of the 
EGMF is computed directly by equating the mag- 
netic energy to a suitable fraction of the turbulent 
energy of the intergalactic gas (see [8] for details). 

Fig. 1 shows the volume fraction in the pgas - B 
plane in our model. The EGMF are correlated with 
the large scale structure: the strongest around clus- 
ters of galaxies and the weakest in voids. In the re- 
gions of galaxy clusters with Pgas/{Pgas) ^ 10^, 
{B) ^ 10~^G. For typical filamentary regions 
with PgasKPgas) ~ 30, {B) - IQ-^G, while in 
void region, {B) - 10" 

CR Sources and Observers 

Assuming that UHECRs originate from AGNs, we 
identify galaxy clusters with kT > 1.0 keV in the 
cosmological simulation as sources, since AGNs 
are likely to form in high density environment. 
There are 20-30 such clusters in the simulated vol- 
ume, corresponding to the source density of = 
2 — 3 X 10~^/i^Mpc"'^ and the mean separation 
of Is ~ 40ft~^Mpc. The field strength in source 
locations ranges 0.1 < Bs{nG) < 10 with (Bg) ~ 



1/iG. UHE protons with a power-law energy spec- 
trum of E-'y for 10 EeV < E < lO^EcV are 
distributed among the sources and then launched 
in random directions. We follow the trajectory of 
individual protons by numerically integrating the 
equation of motion in the model EGMF, consider- 
ing continuous energy loss due to photo-pair and 
photo-pion production [3]. 

An observer is placed at each 10'^ groups of galax- 
ies with 0.05kcV < kT < O.SkcV in order to se- 
lect locations similar to the Local Group. These 
groups are not distributed uniformly, but located 
mostly along filaments, following the matter dis- 
tribution of the large scale structure. Around ob- 
server locations, lO""' < Bobsif^G) < 0.1 with 
(Bobs) ~ lOnG. Once a particle visits an observer 
within a sphere of Robs — 0.5/i"^ Mpc, the arrival 
direction, time delay and energy of the particle are 
registered as an 'observed' event. We let the par- 
ticle continue its journey, visiting several groups 
during its full flight, until its energy is reduced to 
10 EeV. With lO''^ protons injected at sources, about 
10^ events are recorded. 

Results 

Deflection Angle and Time Delay 

Note that the gyroradius of proton is rg — 
lOkpc iE/10^^eY){B/fiG)-^. With our model 
EGMF, the trajectories of ^^10 EeV protons could 
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Figure 2: (a) Distribution of events in the plane of observed particle energy and deflection angle, (b) 
Distribution of events in the plane of observed particle energy and time delay. Protons are injected with a 
power-law spectrum of 7 = 2.7 at sources. 



be severely deflected when they leave the source 
regions with (Bg) ^ 1/iG and when they fly by the 
cluster/group regions with B Si lO^^G, since clus- 
ters/groups have a typical radius of (i?) ~ IMpc. 
According to Fig. 1 the volume filling factor with 
B > IQ-^G is /(> 10"«G) « 0.02. 
In Fig. 2(a) we present the distribution of the de- 
flection angle between the arrival direction of UHE 
protons and the source direction, 9, as a function 
of the observed energy. This distribution depends 
rather sensitively on the power-law index of the in- 
jection spectrum, 7, and the minimim distance be- 
tween sources and observers, D. Here, the results 
of the case with 7 — 2.7 and D = 20Mpc are pre- 
sented. On average the deflection angle decreases 
with the energy, showing a transition from diffuse 
transport regime to rectilinear propagation region 
at - 100 EeV. For E > 40 EeV 30 % of 'ob- 
served' events arrive with 6 < 5° from source di- 
rections, while 60 % with 9 < 5° for E > 100 
EeV. Some super-GZK protons can avoid flying 
through strong field regions around clusters and do 
not get significantly deflected by weak fields in fil- 
amentary structures. 

Since the deflected path length is longer than the 
rectilinear distance, the particle travel time is de- 
layed compared to the light travel time. In Fig. 
2(b), we show the distribution of t^;, the delay of 
the particle travel time compared to the light travel 
time, as a function of the observed energy. The 
rectilinear travel time for the mean separation of 
sources, tree = Is/c ~ lO*'^^ years. On average 



the time delays are longer for lower energy parti- 
cles, with {td) ^ lO^yrs ~ Qtrec around 10 EeV 
and {td) ^ lO^yrs ^ lO'^^^^^ above 100 EeV. 

Energy Spectrum 

Fig. 3 shows the predicted energy spectra of 
UHE proton events for the injection spectrum of 
Nrnj{E) oc E-"' with 7 = 2.0, 2.4 and 2.7 and 
D = 5, 10, and 20 Mpc. We note that J{E)E^ 
in the plot is given in an arbitrary units, since the 
amplitude of the injected spectrum at sources is ar- 
bitrary. In all the cases the presence of GZK su- 
pression above 50 EeV is obvious, so the AGASA 
data cannot be consistent with our results. In case 
of D = b Mpc, there is a significant number of 
particles above 100 EeV that come mainly from 
nearby sources. Below 50 EeV, the spectrum gets 
flatter due to the pileup of the particles that had 
higher initial energies but have lost their energy 
via photo-pion production. The injection spectrum 
with 7 = 2.7 seems to produce the best fit for the 
HiRes-2 data, although the case with 7 = 2.4 can 
be considered as being consistent with the HiRes 
observations. 

Conclusion 

We study the propagation of UHE protons that 
originate from cosmological sources located at 
galaxy clusters, travel through the EGMF corre- 
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Figure 3: Predicted energy spectra of UHE pro- 
tons that are injected with a power-law spectrum 
of Ninj{E) oc E~'^ at sources and propagated 
through a universe with the adopted model EGMF. 
Top panel: the case with the minimum source- 
observer distance D = 20Mpc is shown for 7 = 
2.0 (solid line), 2.4 (dotted) and 2.7 (dot-dashed). 
Bottom panel: the case with 7 = 2.7 is shown for 
D = 5 Mpc (solid line), 10 Mpc (dotted) and 20 
Mpc (dot-dashed). Filled circles are the observa- 
tional data from HiRes-2 [9]. and open circles are 
from HiRes- 1, while stars are from AGASA [7]. 

lated with the large scale structure of the uni- 
verse, and arrive in the regions similar to the Local 
Group. A novel model for the EGMF based on tur- 
bulence dynamo is adopted. Since the volume fill- 
ing factor, /(> 10^ ''G) is less than 5 x lO^'^, about 
60 % of UHE protons above 100 EeV could avoid 
visiting strong field regions and arrive at Earth 
within 5° of their sources, although the deflection 
angle is overall distributed over a wide range. Be- 
low the GZK energy, the propagation is in the dif- 
fusion limit and both the deflection angle and time 
delay are substantial. This indicates that in the 
present scenario, UHE proton astronomy may be 
possible only for > 100 EeV. The predicted en- 
ergy spectrum of UHE protons exhibits the GZK 
supression, and fits the HiRes observations, if the 
injection spectrum has 7 = 2.4 — 2.7. 
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